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Abstract. The study of slope parameter is presented for elastic proton-proton and antiproton-proton scat¬ 
tering with taking into account the resent experimental data at high energies. The expanded logarithmic 
approximations allow the description of the experimental slopes in all available energy range reasonably. 
Accounting for the LHC results leads to the dramatic change of behavior of the quadratic in logarithm 
approximation at high energies and to the closer trends for all fitting functions under study in comparison 
with the analysis at collision energies up to the 200 GeV. The estimations of the asymptotic shrinkage pa¬ 
rameter a'-p are discussed. Predictions for diffraction slope parameter are obtained for some proton-proton 
and antiproton-proton facilities. 

PACS. 13.75.Cs Nucleon-nucleon interactions - 13.85.Dz Elastic scattering 


1 Introduction 

The one of the important quantity for nucleon elastic scat¬ 
tering is the slope parameter B which is defined in accor¬ 
dance with the following equation: 

B{s,t) = dtlndtcr{s,t), (1) 

where dt = djdt. The B is determined experimentally by 
fitting the differential cross section da/dt at some colli¬ 
sion energy a/s. The study of B parameter is important, 
in particular, for reconstruction procedure of full set of 
helicity amplitudes for elastic nucleon-nucleon scattering 
[1]. In the last 20-30 years, high-energy pp colliders have 
extended the maximum pp collision energy from a/s ^ 20 
GeV to ^/s ^ 2 TeV, the LHC facility allows one to ob¬ 
tain pp data up to -/s = 8 TeV so far. As consequence, 
the available collection of pp and pp slope data from litera¬ 
ture has extended. The B{s) for elastic pp and pp reactions 
is under consideration below. Within the classical Regge 
model the Pomeron trajectory, a'p{t), is linear function 
of momentum transfer, i.e. a'p{t) = a-piQ) -I- a'-pt. There¬ 
fore using the dehnition m the following relation can be 
obtained for B{s) at some fixed t: 

B{s) oc 2 q;/ Ine, 

where e = s/sq, sq = 1 GeV^. In general case for Pomeron- 
inspired models the asymptotic shrinkage parameter a'p 
can be written as follows: 2a/(s) = dineB{s,t). 

^ e-mail: VAOkorokov@mephi.ru; Okorokov@bnl.gov 


2 Experimental slope energy dependence 

As in the previous analyses [2113114] dependence B{s) is 
parameterized by the following analytic functions: 


B{s,t) = Bq + 2ai\ne, (2a) 

B{s,t) = Bq + 2ai In e -b 02 ln“^ e, (2b) 

B{s,t) = Bq + 2ai\ne + a 2 S°'^ , (2c) 

i?(s, t) = i?o + 2ai Ine-b 02 In^ e, (2d) 


where free parameters Bq, o^, j = 1 —3 depend on range of 
|t| which is used for approximation. There are the relations 
a'p = oi and a/ = oi -b 02 Ine for parameterizations (l^al) 
and (j2d|) inspired by the Pomeron exchange models. Ap¬ 
pendix shows database of available experimental results 
for slope parameter in elastic pp/pp scattering. The some 
numerical values of the i?(s, t) from this database are used 
in the present study. Experimental values of slope param¬ 
eter collected at initial energies -/s <1.8 TeV were used 
in |3] . Additional experimental results from Tevatron and 
the LHG are from [21] and [MlHmiTTlIH!!] , respectively. The 
full data sample consists of 490 experimental points. The 
number of experimental points is equal 145/138 (137/70) 
for pp/pp scattering at low (intermediate) |t| respectively. 
It should be noted that for intermediate |t| range the ex¬ 
perimental results for B{s) obtained with help of linear 
parametrization for logarithm of differential cross-section, 
\n{da/dt) oc {—B\t\), are discussed below because the 
new experimental data at intermediate |t| [531191] with 
respect to the [4] are obtained for such parametrization 
of In (da/dt) namely. From the exponential parametriza¬ 
tion with index quadratic in t for differential cross-section. 
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In (dcr/dt) cx {—B\t\ one may calculate the local 

slope as following 


Hs,t\^^o,^aeV-=B±C\n\t\, 5,00. (3) 

In accordance with [3] the mean value of |t| (|t|) is calcu¬ 
lated with taking into account the approximation of exper¬ 
imental da/dt distribution; errors of experimental points 
include available clear indicated systematic errors added 
in quadrature to statistical ones. The fitting algorithm is 
described in detail in [3]. As previously the points with 
> 2 are excluded from fit in our algorithm, where [ 3 ] 


f Bj, -B{si;a) \ 
y CTi VxVn-d-f- ) 


( 4 ) 


Here 5^ is the measured value of nuclear slope at Si 
with experimental error cr^, B (s^; a) is the expected value 
from the fitting function with best y^/n.d.f. among (I2al) 
- (IMl) for approximation of all range of available ener¬ 
gies, the parameters aj are from TV-dimensional vector 
a = {ai, ...,ajv}. We consider the estimates of fit param¬ 
eters as the final results if there are no excluded points 
for present data sample. The fraction of excluded points 
is about 2 % for pp as well as for pp elastic scattering for 
low |t| domain. The maximum relative amount of rejected 
points is about 3%/12% for linear hida/dt parametriza- 
tion at intermediate |t| values for pp/pp scattering respec¬ 
tively. 


2.1 Low |f| domain 

The experimental dependence B[s) and corresponding fits 
by ((2a1l - (l 2 dl) are shown in FiglT]for elastic pp scattering. 
Table [1] contains values of ht parameters. As seen the fit¬ 
ting functions (EiD, IMI describe the pp experimental data 
statistically acceptable only for ^/s > 5 GeV. The RHIC 
point at y/s = 200 GeV does not contradict the common 
trend within a large error bars and can’t discriminate the 
approximations under study. In general the LHC results 
iMunmni added in fitting sample result in better agree¬ 
ment of fits (I2a|l “ (I2d|l in comparison with previous study 
[ 3 ]. All fits are very close to each other in energy domain 10 
GeV < y/i < 1 TeV, the quadratic in logarithm function 
(| 2 dp shows some faster increasing of slope and noticeable 
difference from another fit functions in multi-TeV region 
only. It seems the ultra-high energy domain is suitable for 
separation of various parameterizations. Fitting functions 
(Eil), (Ec]) allow us to describe experimental data at all 
energies with reasonable fit quality for pp (Table [T]). The 
functions (1^ - (Eel) agrees very well for y/s > 5 GeV, 
furthermore there is no visible difference between modifi¬ 
cations (I 2 bll and (Eel) in all experimentally available energy 
domain. The function (I2cl) demonstrate some better qual¬ 
ity for fit of full data sample (Table [IJ than function (I2bl) 
in contrast with previous analysis [ 3 ]. The mean value of 
o for excluded points (n^ ) is equal 4.7 for pp data sam¬ 
ple for parametrization (2c). The accounting for LHC data 


leads to some decreasing of values of Bq and increasing of 
ai parameters for all fitting functions (l2al) - (1^ under 
study in comparison with values of corresponding parame¬ 
ters in previous investigation [ 3 ]. This behavior of oi with 
collision energy agrees well with predicted growth of a'^ 
with increasing of [7511M] . As seen from Table [T] the 
third term in both the (I2b|) and the (l2^ gives the main 
contribution at y/s < 5 GeV in the case of elastic proton- 
proton scattering, i.e. describe the sharp changing of slope 
in the low energy domain. Therefore B^^ oc Ine at high 
y/s in accordance with (l2bl) and (E^ . Such asymptotic be¬ 
havior is in qualitative agreement with energy dependence 
of the slope parameter for pp collisions [ 3 ] . The increasing 
of tti exhibits that B^p growth some faster in multi-TeV 
region than one can expect from the trend based on data 
sample at ^/s < 200 GeV. This suggestion is confirmed 
by improvement of fit quality for one fitting function (I2dl) 
for present data sample in comparison with fit qualities 
for experimental points at ^/s < 200 GeV [ 3 ]. Values of 
02 obtained in present study and for fit at ^/s < 200 GeV 
are close within errors for functions (I2bl) and (I2dl) but ac¬ 
counting for the LHC data results in some increasing of 
the absolute value of 03 in fit by ES)- The absolute values 
of 02 and 03 are the same within error bars for function 
(E^ for Table [I] and for fit in energy domain y/s < 200 
GeV g. 



Fig. 1. Energy dependence of the slope parameter for elastic pp 
scattering for low |t| domain. Experimental points from fitted 
samples are indicated as •, unfitted points are indicated as =t=. 
The dot curve is the fit of experimental slope by the function 
(I2all . thick solid - by the ([20, dashed - by the ([20, thin solid 
- by the (I2dl) . The shaded band corresponds to the spread of 
fitting functions for previous analysis [4]. 
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Table 1. Values of parameters for fitting of slope energy dependence in pp elastic scattering for low |f| domain 


Function 



Parameter 




Bo, GeV-2 

ai, GeV^ 

02, GeV^ 

03 

xVn.d.f. 

(I2all 

8.00 ±0.06 

0.309 ±0.010 

- 

- 

234/98 

(f2bl) 

8.09 ± 0.06 

0.305 ± 0.005 

-31.8 ± 1.6 

-4.06 ±0.13 

420/138 

(|2cl) 

7.95 ± 0.06 

0.313 ±0.005 

-240 ± 32 

-2.23 ±0.09 

402/138 

(f2d)) 

8.81 ±0.12 

0.198 ±0.015 

0.013 ±0.002 

- 

174/97 


Table 2. Predictions for slope in nucleon-nucleon elastic scattering at intermediate energies for low |t| domain 


Fitting 



Facility energies ^/s, GeV 


function 


FAIR 

NICA 


3 

5 

6.5 

14.7 20 

25 

(I2al) 

([^ 

(Eel) 

([M]) 

12.57 ±0.07 
12.59 ±0.05 

11.4 ± 0.2 

12.80 ±0.07 
12.83 ±0.06 
11.9 ± 1.2 

11.6 ±0.2 12.3 ±0.2 11.70 ±0.13 

■ 12.93 ±0.07 13.32 ± 0.09 11.72 ±0.08 

; 12.95 ±0.06 13.34 ± 0.07 11.70 ±0.08 

12.0 ±1.3 12.4 ±1.6 11.6 ±0.2 

11.98 ±0.14 

12.00 ±0.09 
11.98 ±0.09 
11.9 ±0.2 

Table 3. Predictions for slope in 

pp elastic scattering at high energies for low |t| domain 

Fitting 



Facility energies ^/s, TeV 


function 

RHIC 


LHC 

FCC/VLHC 



0.5 

14 

28 

42* 100 200 

500 

(I2all 

(120 

(|20 

([20 

15.7 ±0.3 
15.67 ±0.14 
15.73 ±0.14 
15.7 ±0.5 

19.8 ±0.4 
19.7 ± 0.2 

19.9 ±0.2 
21.1 ± 0.9 

20.7 ± 0.4 
20.6 ±0.2 

20.8 ±0.2 
22.4 ± 1.0 

21.2 ± 0.4 22.2 ± 0.5 23.1 ± 0.5 

21.1 ± 0.2 22.1 ± 0.2 23.0 ±0.3 

21.3 ± 0.2 22.4 ± 0.2 23.2 ± 0.3 

23.1 ±1.1 24.8 ± 1.3 26.2 ±1.4 

24.2 ±0.5 

24.1 ±0.3 
24.4 ± 0.3 

28.2 ± 1.6 


‘The ultimate energy upgrade of LHC project [96]. 


The comparison with earlier studies confirms 

the conclusion above that accounting for high energy data 
points leads to increasing of value of parameter and 
growth of slope parameter seems to be faster in TeV-region 
than that at lower energies. Accordingly the value 
for function (iMl) is closer significantly to the a'^ ~ 0.25 
GeV“^ than that for previous analysis in energy range 
•v/s < 200 GeV |4]. On the other hand the values of 
the parameter obtained for fitting function (I2al) is 
some larger than the prediction for a'-p from Pomeron in¬ 
spired model for TeV-energy domain | 94 ] . Furthermore the 
proton-proton results from fit by function (I2dl) allow the 
estimation (s)|= 0.86 ± 0.08 which is almost 
twice larger than the corresponding prediction from [ 94 ] . 

Predictions for B are obtained for some facilities based 
on the fit results shown above for pp and on the results 
from | 3 ] ioT pp. Estimations at low |t| for different interme¬ 
diate energies of the projects FAIR and NICA are shown 
in the Table [5] and for high energy domain are presented 
in the Table |31 As expected the functions (l2bl) and ( 1 ^ 
predicts for FAIR the B values coincide with each other 
within errors. The approximation functions (I2al) and (I2dl) 
can predict for ^/s > 5 GeV only. The Pomeron inspired 
function (l2a1) predicts slope parameter values smaller sig¬ 
nificantly than that for modified fitting functions (I2bl) and 
(E3) at FAIR energies. For pp the predictions with help 


of quadratic in Ine function (I2dl) are equal with estima¬ 
tions based on any another fitting function under study 
within large error bars at y/s < 14.7 GeV. In the case 
of elastic pp scattering the functions (Eil) - (Ecj) predict 
just the equal values of B within error bars at any col¬ 
lision energy ^/s discussed here. The difference between 
estimations from (Pa)) - (23) and from dSl) onsets at the 
final energy of the LHG project = 14 TeV only. Our 
prediction with (I2dl) function for RHIC energy is equal 
with early prediction for close energy based only on slope 
data in the region 5 < i/s < 62 GeV [37] within errors. 
But (l2dl) underestimates the B values in ultra-high en¬ 
ergy domain ^/s > 40 TeV in comparison with results 
based on the approach without odderons m- It should 
be emphasized that in contrast with previous analysis [3] 
the present fits by functions (l2al) - (iMll of data sample 
included of the LHG results predict the similar increas¬ 
ing of B with energy as most of phenomenological mod¬ 
els (98). The B value predicted for the LHG at y/s = 14 
TeV by (l^al) only is close with errors to the predictions 
from [MlfTHn) . in particular, with estimation from model 
with hadronic amplitude corresponding to the exchange 
of two pomerons. Prediction of phenomenological model 
with hadronic amplitude corresponding to the exchange 
of three pomerons [100] at y/s = 14 TeV coincides with 
estimation of B within error bars from fit function (I2dl) 
with fastest growth of B with ^/s in multi-TeV region. But 
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most of estimations of B at y/s = 14 TeV from Table [3] 
agree well within errors with model prediction from m- 
However the model estimates at y/s = 14 TeV described 
above were obtained for B {t = 0) and the t-dependence 
of slope shows the slight decreasing of B for the model 
with three-pomeron exchange |100) and faster decreasing 
of B for the model from m at growth of momentum 
transfer up to |t| ~ 0.1 GeV^. Therefore one can expect 
that the model with hadronic amplitude corresponding to 
the exchange of three pomerons m will be in the bet¬ 
ter agreement with values of B from Table [3] predicted for 
finite (non-zero) low |t| values. But the model from |102) 
overestimates the B value at y/s = 14 TeV in comparison 
with corresponding predictions from fitting functions (I2al) 
- (Edl) despite of the sharp decreasing of B at growth of 
momentum transfer up to |t| Ri 0.1 GeV^ in this model. 
As suggested sometimes the saturation regime, Black Disk 
Limit (BDL), maybe reach at the LHG. The one of the 
models in which such effects appear, namely, Dubna Dy¬ 
namical Model (DDM) predicts the slope B (t = 0) Ri 23.5 
GeV“^ at ^/s = 14 TeV [103] which is larger noticeably 
than the predictions from Table [3] at the same ^/s. There¬ 
fore the saturation regime will not be reached, at least, at 
the LHG energy = 14 TeV as suggested, for example, 
in the model from [104] or simple saturation can not be 
enough in order to describe the LHG data at quantitative 
level. 


2.2 Intermediate \t\ domain 

As indicated in [3] the situation is more complicated for 
intermediate |t| domain. Figure [2] shows the experimental 
data and corresponding fits for energy dependence of slope 
parameter at intermediate |t| for exponential approxima¬ 
tion with index linear in |t| of differential cross-sections in 
elastic pp and pp scattering. The parameter values for fit¬ 
ting are indicated in Table [3] for various interaction types. 

Usually the fit qualities are poorer for intermediate |t| 
values than that for low |t| range in pp elastic collisions 
for linear parametrization of ln{d(T/dt). The fitting func¬ 
tions (I2al) and (I2dp agree with experimental points qual¬ 
itatively for -^/s > 5 GeV only. Furthermore in the first 
case there is significant discrepancy between experimen¬ 
tal point and fit curve at the LHC energy also (Figl3^). 
The “expanded” functions (l2bl) . (I3cl) approximate experi¬ 
mental data at all energies reasonably with close fit qual¬ 
ities (Table |3|), but these functions show a slow growth 
of slope parameter with energy increasing at ■\/s > 100 
GeV (Fig|2^). It should be stressed that the experimental 
point at the LHC energy leads to the dramatic change of 
behavior of the fitting function (I2dl) in comparison with 
previous analysis [3]. At present the fitting function (I2dl) 
predicts increasing of the nuclear slope in high energy do¬ 
main as well as all other fitting functions under study. 
Such behavior is opposite to the result of fit by function 
((Ml) of experimental data sample at y/s < 200 GeV [3]. 
As shown in previous analysis [3] the the pp experimental 
data admit the approximation by (I3al) . (IMl) for all en¬ 
ergy range but not only for ^/s > 5 GeV (Figl^lo). Thus 


the parameter values are shown in Table [3] for approxi¬ 
mation by (l2al) . (I2dl) of all available experimental data. 
The fit curves show (very) close behaviors for both the 
present and previous analyses in the case of the Figl23. 
The pp data disagreement with Pomeron inspired fitting 
function (I3al) very significantly (Figl3)D). Functions (IMl) 
and (l2cll show a close behavior at all energies for pp data 
from linear parametrization of In da/dt. These fitting func¬ 
tions have a better fit quality than (IMl) but fits by func¬ 
tions (|2bll . (|2cll are still statistically unacceptable. The 

= 2.9 for excluded pp data with (I3cl) function and 

= 18.3 for points excluded from pp fitted data sam¬ 
ple for (I2bll fitting function. As seen from Figl^Jo behavior 
of the energy dependence of slope parameter for elastic 
pp scattering is close to the quadratic in logarithm func¬ 
tion BPP oc In^ £ at high ^/s as well as for elastic pp col¬ 
lisions. But in the last case it is difficultly to make the 
unambiguous conclusion because there is only the LHC 
point in the asymptotic region. The estimations of asymp¬ 
totic shrinkage parameter are following for fitting func¬ 
tions (l2ail and (IMll respectively: 2a'-p = 0.422 ±0.016 and 
2a)),(s)|^^.j,rpgy = 1.66±0.19 for pp; 2a'j, = 0.272±0.008 
and 2ap(s)|^^^ gexeV ~ ^ PP- values 

of a'-p are equal within errors at ^/s = 7 TeV and ^/s = 8 
TeV for domain of |f| under consideration in the case of 
the function (I2dll for pp interaction. As seen from compar¬ 
ison with results from Table [1] for elastic pp collisions and 
from [3] for pp scattering the fitting function (I2al) is char¬ 
acterized by smaller asymptotic shrinkage parameter for 
intermediate |t| domain than that at small |t| values while 
the situation for quadratic in logarithm function (I2dl) is 
opposite: there is significant growth of a'p for transition 
from small \t\ to intermediate \t\ values. 

We have obtained predictions for nuclear slope param¬ 
eters B and b for some facilities and intermediate |<|. The 
B values are calculated on the base of fit results shown 
above. Estimations for b are obtained on the base of results 
from [4] with all functions under consideration with excep¬ 
tion of the (I2bl) . In the last case the following results are 
obtained within framework of the present analysis for 6(s): 
Bo = 9±2 (50.8T0.9) GeV-^ oi = 0.19T0.12 (0.92T0.06) 
GeV-2, 02 = -10 ± 2 (-34.1 ± 0.8) GeV-^, og = -1.9 ± 
0.5 (0.222T0.016), xVn.d.f. = 41.4/33 (8.7/15) for pp (pp) 
data samples, respectively. Slope values are shown in the 
Table [5] for different energies of FAIR, NICA, and in the 
Table [6] for RHIG, LHG and FCC/VLHC. According to 
the ht range function (I3al) can predicts the B value for pp 
scattering at all energies under study not in > 5 GeV 
domain only. As expected the functions (IM)) - (IMll pre¬ 
dicted very close values for slope parameter B for FAIR. 
These fitting functions, especially ((20 and (I3cll . predict 
the close values for nuclear slope B in NIGA energy do¬ 
main too. Functions (I3ai) - ((3cl) predict smaller values for 
B in high-energy pp collisions than (I2dl) approximation 
especially for FGG/VLHG energy domain. Perhaps, the 
future more precise RHIG results will be useful for dis¬ 
crimination of fitting functions under study for interme¬ 
diate |t| values. In the contrast with previous analysis [3], 
here the function dMl with obtained parameters predicts 
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Fig. 2. Energy dependence of B in proton-proton (a) and antiproton-proton (b) elastic scattering for exponential parametriza- 
tion with index linear in |t| of differential cross section. Experimental points from fitted samples are indicated as close circles / 
squares for pp/pp, unfitted points are indicated as The dot curve is the fit of experimental slope by the function (l2aD . thick 
solid - by the (123, dashed - by the (123, thin solid - by the ([ 23 . The shaded band (a) corresponds to the spread of fitting 
functions for previous analysis [3]- 


Table 4. Values of parameters for fitting of energy dependence of slope at intermediate |t| 
Function Parameter 

Bo, GeV~^ ai, GeV~^ 02, GeV~^ 03_xVn.d.f. 


proton-proton scattering, experimental data for da/dt = Aexp(—i3|f|) 


(I2all 

7.59 ± 0.11 

0.211 ±0.008 

- 

- 

322/35 

(HB 

8.39 ±0.17 

0.163 ± 0.011 

-25.2 ±1.4 

-3.01 ±0.13 

493/61 

(123 

7.94 ± 0.11 

0.19 ±0.09 

-90 ±8 

-1.69 ±0.06 

458/61 

(123 

9.9 ± 0.2 

-0.16 ±0.03 

0.056 ± 0.005 

- 

187/34 


antiproton-proton scattering, experimental data for da/dt = j4exp(—i3|t|) 


([2a 

11.16 ±0.06 

0.136 ±0.004 

- 

- 

1168/42 

(123 

14.34 ±0.11 

-0.304 ±0.014 

0.0042 ± 0.0004 

2.92 ±0.04 

186/40 

(I 23 

7.7 ±0.2 

-0.735 ±0.019 

6.9 ±0.2 

0.1000 ±0.0013 

188/40 

(123 

15.48 ±0.15 

-0.60 ±0.02 

0.084 ± 0.003 

- 

197/41 


Table 5. Predictions for slopes in nucleon-nucleon elastic scattering at intermediate energies for intermediate |t| domain 


Fitting 



Facility energies ^/s, GeV 



function 


FAIR 


NICA 


3 

5 

6.5 

14.7 

20 

25 

B-parameter 

(2a) 

11.76 ±0.06 

12.04 ±0.07 

12.18 ±0.07 

12.62 ±0.07 

10.12 ±0.15 

10.31 ±0.15 

(2b) 

12.88 ±0.13 

12.51 ± 0.14 

12.26 ±0.15 

11.6 ±0.2 

10.2 ±0.2 

10.4 ± 0.2 

(2c) 

13.1 ±0.3 

12.5 ±0.4 

12.2 ±0.4 

11.6 ±0.5 

10.2 ± 1.1 

10.4 ± 1.2 

(2d) 

13.27 ±0.18 

12.5 ±0.2 

12.2 ±0.2 

11.5 ±0.3 

10.0 ±0.5 

10.1 ±0.5 

6-parameter 

(2a) 

10.0 ±0.2 

10.5 ±0.2 

10.7 ±0.3 

11.4 ± 0.3 

10.5 ±0.3 

10.8 ±0.4 

(2b) 

14.3 ± 1.4 

12.6 ± 1.7 

12.0 ± 1.8 

11 ± 2 

10 ±2 

11 ± 3 

(2c) 

14 ±3 

13 ±3 

12 ±3 

11 ± 3 

10.5 ± 1.1 

10.7 ± 1.2 

(2d) 

12.6 ± 1.1 

12.0 ± 1.2 

11.8 ± 1.3 

11.3 ± 1.7 

10 ±2 

11 ± 3 
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Table 6. Predictions for slopes in pp elastic scattering at high energies for intermediate |t| domain 


Fitting 

function 



Facility energies y/a, TeV 



RHIC 


LHC 

FCC/VLHC 



0.5 

14 

28 42* 

100 

200 

500 

B-parameter 

(2a) 

12.8 ±0.2 

15.6 ±0.3 

16.2 ±0.3 16.6 ± 0.4 

17.3 ±0.4 

17.9 ±0.4 

18.7 ± 0.4 

(2b) 

12.4 ±0.3 

14.6 ±0.4 

15.1 ±0.5 15.3 ±0.5 

15.9 ±0.5 

16.3 ±0.6 

16.9 ±0.6 

(2c) 

13 ±2 

15 ±3 

16 ±4 16 ±4 

17 ±4 

17 ±4 

18 ±5 

(2d) 

14.5 ± 1.1 

24 ± 2 

27 ± 2 28 ± 3 

32 ±3 

35 ±3 

40 ±4 

6 -par ameter 

(2a) 

14.6 ±0.6 

18.8 ±0.9 

19.7 ± 1.0 20.2 ±1.0 

21.3 ± 1.1 

22.2 ± 1.1 

23.3 ± 1.2 

(2b) 

13 ±4 

16 ±5 

16 ±5 17 ±5 

17 ±6 

18 ±6 

18 ±6 

(2c) 

13.7 ± 1.9 

17 ±3 

18 ± 3 18 ± 3 

19 ±3 

20 ±3 

21 ±4 

(2d) 

12 ±5 

9±9 

8 ±10 7 ±11 

6 ± 12 

4± 13 

2± 14 


*The ultimate energy upgrade of LHC project Oil- 


fast growth of B values at energies of future experiments. 
This behavior of estimations calculated for functions (I2al) 
- (I2dl) contradicts with earlier predictions from some phe¬ 
nomenological models. It should be emphasized that var¬ 
ious phenomenological models predict a very sharp de¬ 
creasing of nuclear slope in the range \t\ ^ 0.3 — 0.5 GeV^ 
at LHC energy y/s = 14 TeV [SH]. Just the positive B 
value predicted for LHC at y/i = 14 TeV by (IMll is most 
close to the some model expectations fTIMfinTI . Taking 
into account predictions in Table [H based on the fitting 
functions (^3 - (l2d)) for low |t| one can suggest that the 
model with hadronic amplitude corresponding to the ex¬ 
change of three pomerons piirnn] describes the nuclear 
slope some closer to the experimentally inspired values at 
LHC energy both at low and intermediate |t| than other 
models. The situation with predictions for b at interme¬ 
diate energies is similar to that for B: functions (I2bl) - 
(I2dl) predict close values for b within large errors for both 
the FAIR and the NICA. Furthermore all fitting functions 
predict the close values of b at y/s > 5 GeV. The value of b 
parameter obtained from function (I2all differs significantly 
from estimations with other fitting functions for pp elas¬ 
tic scattering at low energies. Therefore the b parameter 
seems more perspective for distinguishing of Pomeron in¬ 
spired function (I2al) from “expanded” parameterizations 
dSi) - © at ^/s ^ 3 GeV than B. One can see the 
functions (l2bll and (l2^ predict very close values of b in 
high-energy domain. Function (I2dl) shows a decreasing of 
b at high energies in the contrast with B parameter. In 
general estimations obtained with help of all fitting func¬ 
tions are agree within errors up to the y/s = 28 TeV. But 
the large errors for function (j2dll do not allow the unam¬ 
biguous physics conclusion especially at the LHC energies 
and above. 


2.3 AB and NN data analysis 

In accordance with rules from [3] the difference of slopes 
(AB) for antiproton-proton and proton-proton elastic scat¬ 
tering is calculated for each function (Eal) - (l2d)) under 
study with parameters from corresponded pp and pp fits: 


ABi{s) = Bf^{s) — Bf^{s), i = (2a) - (2d). The energy 
dependence of AB is shown at Figl3^ and Figl3}D for low 
and intermediate |<| respectively. One can see that the dif¬ 
ference of slopes decreasing with increasing of energy for 
low \t\ domain (Figl3^) as well as in the previous analy¬ 
sis [1]. The fitting functions (I2bl) . (1^ demonstrate much 
faster decreasing of AB with increasing of y/s than that 
the functions (I2a1l and (iMll . At present the proton-proton 
experimental data at highest available energy 8 TeV don’t 
contradict with fast (square of logarithm of energy) in¬ 
creasing of slope at high energies in general case. Such 
behavior could be agreed with the asymptotic growth of 
total cross section. Furthermore in contrast with the pre¬ 
vious analysis [3], here the quadratic in Ine function (IMl) 
leads to much smaller difference AB for pp and pp scat¬ 
tering in high energy domain for both low (Fig|3^) and 
intermediate (FigOb) values of |t|. The Pomeron inspired 
function (l^all only predicts the decreasing of AB with en¬ 
ergy growth at intermediate t\ (FiglSja) for any values of 
^/s. The parameterizations (I2bl) - (I2dl) predict the decreas¬ 
ing of difference of slopes at low and intermediate energies 
and fast increasing of AB at high energies for intermediate 
|t| domain (FigjSjD). As expected the most slow changing 
of AB is predicted by Pomeron inspired function (1^ at 
asymptotic energies. All fitting functions with experimen¬ 
tally inspired parameters don’t predict the constant zero 
values of AB at high energies. But it should be empha¬ 
sized that only separate fits were made for experimental 
data for pp and pp elastic reactions above. These results 
indicate on the importance of investigations at ultra-high 
energies both pp and pp elastic scattering for study of 
many fundamental questions and predictions related to 
the general asymptotic properties of hadronic physics. 

Also we have analyzed general data samples for pp 
and pp elastic scattering. Slope parameters {B and b) 
shows a different energy dependence at y/s < 5 GeV in 
proton-proton and antiproton-proton elastic reactions in 
any |t| domains under study. Thus slopes for nucleon- 
nucleon data are investigated only for y/s > 5 GeV be¬ 
low. We have included in fitted samples only pp and pp 
points which have been included in corresponding final 
data samples at separate study pp and pp elastic reactions 
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Fig. 3. The energy dependence of the difference of elastic slopes for antiproton-proton and proton-proton scattering in low 
|t| domain (a) and in intermediate |f| range for exponential fit with index linear in |t| of differential cross-section (b). The 
correspondence of curves to the fit functions is the same as well as in Figure 1. 




\/s, GeV \/s, GeV 

Fig. 4. Energy dependence of the elastic slope parameter in nucleon-nucleon scattering for low (a) and intermediate (b) domain 
of absolute values of momentum transfer. Experimental points from fitted samples are indicated as • for pp and as ■ for pp. 
Fits are shown for = 20 GeV. The dot curve is the fit of experimental slope by the function (I2all . thick solid - by the 

(120, dashed - by the (pell. thin solid - by the ( 1231 ). The shaded band corresponds to the spread of fitting functions for previous 
analysis H]. 
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above. We did not exclude any points from NN data sam¬ 
ple, we change only the low energy boundary for fitted 
domain. Experimental data for slope in nucleon-nucleon 
elastic scattering against collision energy are shown in 
FigSh- at low \t\ and in FigHJo for intermediate |t| together 
with fits by functions (Eal) - (iMl) . We have fitted the gen¬ 
eral nucleon-nucleon data sample at range of low energy 
boundary Sjjjjjj = 25,100, 225 and 400 GeV^. The ht pa¬ 
rameter values are indicated in Table [7] on the first line for 
low boundary of the htted energy domain Smin = 25 GeV^ 
and on the second line - for Sjjjjn = 400 GeV^. The ht 
quality improves for most parameterizations under consid¬ 
eration at increasing of Smini ttius htting functions (l2al) - 
dMI are shown at Fig|3]for = 400 GeV^. As seen from 
FigSh- there is no experimental data for pp between y/s = 
5 GeV and ^/s = 10 GeV at low \t\. This energy domain 
will available for further FAIR facility. One need to empha¬ 
size the ht quality is some poorer (y^/n.d.f. ~ 2.3 — 2.9) 
at ^/s > 10 GeV than that for ^/s > 5 GeV for functions 
(1^ and (I7B . For all cases of Sjjjjjj indicated above the 
value of the oi parameter obtained with the function (I2al) 
agrees qualitatively with the prediction within the frame¬ 
work of Pomeron model, but the value of the asymptotic 
shrinkage parameter {2a'-p = 0.662 ± 0.010) obtained at 
Smin = 400 GeV^ is some larger than the prediction for 
a'-p from Pomeron inspired model for TeV-energy domain 
[M] . Also results from ht by function (I2dl) with accept¬ 
able quality at Sj^in = 400 GeV^ allow the estimation 
2Q!p(s)|^^g.pgy = 0.74 ± 0.12 which is some larger than 
the corresponding prediction from |94) . Furthermore the 
estimations for 2Q!p(s)|^^grpgy do not depend on Smin 
within error bars. At low |t| all functions (I2al) - (I2dp are 
close to each other at energies up to y/i ^ 1 TeV at least 
and shows quasi-linear behavior for parameter values ob¬ 
tained by hts with s^in = 25 GeV^ and Sjjiin = 400 GeV^. 
This observation conhrms the suggestion that oc In e 

at high ^/s at low |t| values. As seen from comparison 
between the present hts and the spread of previous ht 
functions (shaded band) there is a dramatic change of be¬ 
havior of the htting function (I2dl) in comparison with pre¬ 
vious analysis [3] due to experimental points at the LHG 
energies. At present the htting function (l2dl) predicts in¬ 
creasing of the nuclear slope in high energy domain as well 
as all other htting functions under study. Such behavior is 
opposite to the result of ht by function dSI of experimen¬ 
tal data sample at y/s <1.8 TeV [4]. We have analyzed 
the nucleon-nucleon data for slope parameter B at inter¬ 
mediate |<| values for exponential parametrization with 
index linear in |t| of da/dt (FigHJa). As seen experimental 
pp and pp data for B differ signihcantly up to y/s ~ 10 
GeV at least that results in unacceptable ht qualities for 
all functions under study (y^/n.d.f. ~ 28.9 for best ht by 
quadratic in logarithm function). The Pomeron inspired 
function (I2all contradicts with experimental data at any 
Smin • We would like to emphasize that the best ht quality 
for (1^ is obtained at Smin = 100 GeV^ (y^/n.d.f. ~ 9.45) 
but it is statistically unacceptable too. Functions (l2bll - 
dsai) agree with experimental dependence B{s) reasonably 
and have a close ht qualities. Furthermore the functions 


120 and (I2c|l demonstrate very close behaviors in all range 
of y/s under consideration. Best ht is ” expanded” function 
(1^ in difference with previous analysis [3]. As seen from 
Figl3}3 the slope parameter B increases much faster that 
Ine at intermediate |t| values. This result conhrms the sug¬ 
gestion made on the base of the Figl3] above and allows 
us to derive the asymptotic relation B^^ oc In^ e, fc > 2 
at high y/s for intermediate |t| domain. Fig|4)D shows that 
the experimental point at the LHG energy leads to faster 
increasing of most of htting functions in multi-TeV region 
in comparison with previous analysis [ 3 ] . Comparison be¬ 
tween Fig|3h' and Figl3ti demonstrates that the domain 
of intermediate |t| is preferable with respect to the low 
|t| range for discrimination of various phenomenological 
parameterizations for B{s), in particular, the Pomeron in¬ 
spired function (1^ from ’’expanded” ones. At Sjjjjjj = 400 
GeV^ and intermediate |t| values the following estimations 
of asymptotic shrinkage parameter are obtained for htting 
functions (l2a|l and (l2^ respectively: 2a'-p = 0.706 ±0.012 
and 2ap(s)|^^7rpgy = 1.9± 0.2. As well as for low \t\ do¬ 
main the values of a/, are equal within errors at y/s = 7 
TeV and y/s = 8 TeV in the case of the function (iMll : the 
estimations for “^ci'-pis/^^YTeV depend on Sjjjjjj 

within error bars. The growth of a/, is observed with in¬ 
creasing of \t\ (|t|) values, especially, for htting function 

(I2d|. 

One can conclude the slope parameters for pp and pp 
elastic scattering show universal behavior at y/s > 20 GeV 
and “expanded” functions represent the energy depen¬ 
dence for both the low and the intermediate |t| ranges for 
this energy domain. Thus quantitative analysis of slopes 
at different |t| allows us to get the following estimation 
of low energy boundary: y/s ~ 20 GeV for universal¬ 
ity of elastic nucleon-nucleon scattering. This estimation 
agrees qualitatively with both the results for differential 
cross-sections of pp and pp elastic reactions based on the 
crossing-symmetry and derivative relations [T] and the re¬ 
sults for global scattering parameters |104] . 


3 Conclusions 

The present status of diffraction slope parameter for elas¬ 
tic pp and pp scattering is analyzed over the full energy 
domain as well as predictions for some facilities. The ’’ex¬ 
panded” parameterizations allow us to describe experi¬ 
mental B at all available energies in low |t| domain for pp 
as well as for intermediate |t| values for pp and pp quite 
reasonably. The similar situation is observed for hts of 
data samples joined for elastic NN scattering. Therefore 
“expanded” functions can be used as a reliable hts for wide 
range of momentum transfer at all energies. The new LHC 
data lead to dramatic change of behavior of quadratic in 
logarithm function and usually to better agreement be¬ 
tween various htting function in comparison with the anal¬ 
ysis of pp data at y/s < 200 GeV. At low values of |t| the 
’’standard” approximation differs from expanded ones and 
experimental data mostly in the low energy region; at in¬ 
termediate |t| this function is unacceptable for htting of 
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Table 7. Values of parameters for fitting of slope energy dependence in NN elastic scattering 


Function 



Parameter 



Bo, GeV-2 

ai, GeV"^ 

02 , GeV-2 

03 

xVn.d.f. 

low |t| domain 

(|2a|l 

8.10 ±0.05 

0.301 ± 0.004 

- 

- 

332/127 


7.55 ±0.05 

0.331 ± 0.005 

- 

- 

103/63 

(EB 

9.4 ± 0.2 

6.1 ± 1.1 

-12 ±2 

-3.5 ± 1.5 

281/125 


8.1 ±0.7 

2.532 ± 0.006 

-4.6 ±0.2 

0.988 ±0.015 

103/61 

(EB 

8.15 ±0.09 

0.208 ±0.010 

0.50 ±0.09 

0.118 ±0.005 

281/125 


7.8 ±0.2 

0.323 ± 0.009 

-7.0 ± 1.7 

-0.65 ±0.15 

102/61 

EB 

8.81 ±0.12 

0.204 ±0.015 

0.011 ±0.002 

- 

286/126 


8.0 ±0.3 

0.28 ± 0.04 

0.005 ± 0.003 

- 

101/62 

intermediate |t| domain 

(|2a|l 

8.25 ± 0.06 

0.214 ±0.004 

- 

- 

3135/57 


5.17 ±0.11 

0.369 ± 0.006 

- 

- 

353/36 

EB 

12.12 ±0.10 

-0.217 ±0.11 

0.0097 ± 0.0009 

2.60 ± 0.04 

1615/55 


10.7 ±0.2 

-0.019 ±0.011 

(1.04 ±0.17) X 10"'‘ 

4.08 ±0.08 

107/34 

(1^ 

4.43 ±0.17 

-0.602 ±0.013 

7.72 ±0.16 

0.0911 ±0.0009 

1613/55 


10.5 ±0.2 

-0.165 ±0.019 

0.66 ±0.06 

0.187 ±0.004 

106/34 

EU 

13.07 ±0.14 

-0.438 ±0.017 

0.0743 ± 0.0019 

- 

1617/56 


14.9 ±0.6 

-0.61 ±0.06 

0.089 ± 0.006 

- 

118/35 


experimental NN data at all. The intermediate |t| range 
is preferable with respect to the low |t| values for discrim¬ 
ination of various phenomenological parameterizations for 
B{s) dependence. Based on the nuclear slope the low en¬ 
ergy boundary for universality of elastic nucleon-nucleon 
scattering is estimated as ^/s ~ 20 GeV for both the low 
and the intermediate |t| values that is in agreement with 
the hypothesis of a universal shrinkage of the hadronic 
diffraction cone at high energies. The difference of slopes 
for antiproton-proton and proton-proton elastic scatter¬ 
ing {AB) shows the opposite behaviors at high energies for 
low and intermediate |t| domains (decreasing / increasing, 
respectively) for all fitting functions with the exception 
of Pomeron inspired one. All underlying fitting functions 
with experimentally inspired values of parameters don’t 
predict the zero value for AB at both the low and the 
intermediate |t| ranges at high energies. Slop analysis of 
joined NN data samples allows us to estimate the asymp¬ 
totic shrinkage parameter parameter for various domains 
of \t\. The estimation of the a'-p obtained with quadratic 
in logarithm function for NN data at ^/s = 8 TeV for 
low |t| is noticeably larger than the expectation from the 
Pomeron theory. But the growth of a'p with increasing of 
y/s is observed from the comparison of the fit results from 
present study and our earlier analysis for y/s <1.8 TeV. 
The such behavior of a'p agrees with Pomeron inspired 
model. The growth of a'p is observed with increasing of 
momentum transfer. Based on the fit results the predic¬ 
tions for slope parameters B and b are obtained for elastic 
pp and pp scattering in energy domains of some facilities. 
It seems the phenomenological model with hadronic am¬ 
plitude corresponding to the exchange of three pomerons 
describes the B some closer to the experimental fit in¬ 
spired values at the LHC energy both at low and interme¬ 
diate |t| than other models. 


Appendix 

Experimental database for slope in nucleon- 
nucleon elastic scattering 

In this Appendix, the numerical values used in the 
Sec. [5] are shown for diffractive slope parameter. The da¬ 
tabase includes the experimental values of the slope ob¬ 
tained with exponential approximations of da/dt shown 
above. In the Tables 151 - ITT] the statistical errors are shown 
by first and available systematic uncertainties are demon¬ 
strated by second. The compilation of some early experi¬ 
mental measurements of the slope in pp/pp elastic scatter¬ 
ing at low energies was made in [Sj. This one reference is 
shown and detail list of corresponding original papers can 
be found in [^ . In the case of pp elastic scattering in low 
|t| domain (Table [ini) the average y/s are demonstrated 
for ranges of initial momentum (piab) 1-23 - 1.33 GeV/c 
and 1.45 - 1.65 GeV/c [S]. As for database of collective 
parameters |93) the results from latest paper of the exper¬ 
iment are included in Tables below at certain y/s 

and range of |t|. Measurements of particular experiment 
can be separated among Tables [8] - [TT] in accordance with 
the limits for low / intermediate |t| domain dehned above. 
But here there is no averaging over results of various ex¬ 
periments at equal collision energy in difference with |93) . 
Furthermore the averaging over results in different ranges 
of |t| for particular experiment at the same y/s is absent 
too within the low / intermediate |t| domain. Thus the 
measurements of particular experiment at fixed y/s and 
in various ranges of |t| are shown separately in the cor¬ 
responding Table below if these measurements are in the 
chosen |t| domain (see, for example, P^r74l[75] l and agree 
with general trend, at least, qualitatively. Therefore the 
database includes the differential measurements over both 
the y/s and the |t| in the cases of some experiments. 
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V.A. Okorokov: Energy dependence of slope parameter in elastic nucleon-nucleon scattering 


Table 8 . Experimental Bpp for low |t| domain 


yi, GeV 

Bpp, GeV-2 

Ref. 

Vs, GeV 

Bpp, GeV-2 

Ref. 

Vs, GeV 

Bpp, GeV-2 

Ref. 

1.962 

1.60 ±0.60 

[5] 

6.977 

10.70 ± 0.40 

m 

19.42 

10.80 ±0.50 

HU 

1.974 

1.43 ±0.41 

m 

7.167 

10.20 ±0.20 

m 

19.42 

11.50 ±0.40 

HU 

1.993 

0.40 ± 0.50 

[5] 

7.311 

10.52 ±0.12 ±0.30 

m 

19.46 

11.93 ±0.05 ±0.32 

HU 

2.000 

0.38 ±0.19 

m 

7.620 

10.61 ±0.27 

m 

19.66 

12.10 ± 1.20 

HU 

2.002 

0.60 ±0.20 

[5] 

7.664 

10.49 ±0.12 ±0.30 

m 

21.22 

11.61 ±0.19 ±0.20 

cu 

2.069 

0.99 ±0.38 

m 

8.019 

10.69 ±0.12 ±0.30 


21.49 

11.57 ±0.03 

HU 

2.069 

1.46 ±0.47 


8.345 

10.57 ±0.11 ±0.30 


21.70 

12.20 ±0.30 

HU 

2.083 

0.64 ±0.29 


8.550 

10.68 ±0.09 ±0.30 


22.55 

11.69 ±0.10 ±0.20 

CU 

2.125 

2.63 ±0.37 


8.832 

10.82 ±0.11 ±0.30 


23.50 

11.80 ±0.30 

HU 

2.177 

3.56 ± 1.09 

m 

9.030* 

13.50 ±0.74 

m 

23.60 

11.80 ±0.30 

HU 

2.218 

4.42 ±0.51 


9.303 

10.90 ±0.09 ±0.30 

m 

23.76 

12.40 ± 0.30 

HU 

2.239 

5.83 ±0.39 


9.777 

10.50 ± 0.40 

HU 

24.23 

11.90 ±0.28 ±0.20 

CU 

2.259 

4.88 ±0.45 


9.777 

11.25 ±0.28 

m 

24.30 

11.40 ±0.50 ±0.07 

HU 

2.279 

5.53 ±0.46 


9.778 

10.70 ±0.18 ±0.20 

m 

25.59 

11.96 ±0.15 ±0.20 

CU 

2.300 

6.49 ±0.60 


9.778 

10.00 ±0.20 ±0.15 

HU 

26.42 

11.87 ±0.15 ±0.20 


2.319 

6.24 ±0.37 


9.837 

10.84 ±0.11 ±0.30 

HU 

27.29 

11.77 ±0.10 ±0.20 


2.465 

5.70 ± 0.60 

[5] 

9.977 

11.00 ±0.12 ±0.30 


30.48 

12.00 ±0.20 

m 

2.465 

5.97 ±0.15 

[8] 

9.987* 

12.70 ±0.70 

HU 

30.60 

12.00 ±0.20 

HU 

2.696 

7.80 ± 0.50 

m 

10.19 

11.12 ±0.13 ±0.30 

HU 

30.60 

12.20 ±0.30 

HU 

2.702 

7.60 ± 0.43 

[9] 

10.43 

11.11 ±0.10 ±0.30 


30.70 

12.90 ± 1.10 

m 

2.725 

7.80 ± 0.37 

m 

10.52 

10.83 ±0.07 ±0.20 

m 

30.80 

12.30 ±0.30 

m 

2.735 

6.50 ±0.60 

[5] 

10.71 

11.05 ±0.08 ±0.30 

HU 

30.81 

11.87 ±0.28 

m 

2.987 

8.33 ±0.37 


11.00 

11.50 ±0.11 ±0.30 

m 

30.82 

13.00 ± 0.70 

HU 

3.266 

7.69 ±0.37 

m 

11.17 

11.50 ±0.11 ±0.30 

HU 

31.00 

11.70 ±0.62 

m 

3.307 

7.80 ± 0.44 

[9] 

11.22 

11.24 ±0.11 ±0.30 


44.70 

12.80 ±0.30 

m 

3.321 

7.70 ± 0.50 

m 

11.47 

11.46 ±0.09 ±0.30 


44.87 

12.87 ±0.20 

HU 

3.843 

9.60 ±0.30 


11.53 

11.48 ±0.15 ±0.30 


44.90 

13.30 ±0.30 

HU 

3.851 

9.14 ±0.35 

[S] 

11.54 

11.21 ±0.40 

m 

45.00 

12.60 ± 0.40 

HU 

3.874* 

8.20 ±0.20 


11.54 

10.50 ±0.30 ±0.16 

HU 

45.00 

12.80 ±0.30 

HU 

3.875 

9.40 ± 0.60 

m 

11.54 

10.87 ±0.14 ±0.34 

HU 

45.06 

12.90 ± 0.40 

HU 

3.995 

9.53 ±0.23 

[S] 

11.94 

10.80 ±0.50 

HU 

52.80 

13.10 ±0.20 

HU 

4.305 

9.52 ±0.26 

[5] 

12.17 

10.84 ±0.20 ±0.20 

m 

52.80 

13.10 ±0.20 

HU 

4.307 

9.40 ± 0.30 

m 

13.76 

11.60 ±0.50 

HU 

52.80 

12.70 ±0.20 

HU 

4.329 

8.72 ±0.38 ±0.20 

m 

13.76 

11.20 ±0.60 

HU 

52.80 

13.09 ±0.37 ±0.21 

HU 

4.478 

9.50 ±0.22 

[S] 

13.76 

11.13 ±0.14 ±0.33 

HU 

52.80 

12.87 ±0.14 

HU 

4.540 

9.20 ±0.40 

m 

13.76 

12.00 ±0.30 

HU 

52.90 

13.10 ±0.30 

HU 

4.562 

10.40 ± 0.40 

m 

13.90 

11.24 ±0.13 ±0.20 

m 

52.99 

12.40 ± 0.30 

HU 

4.621 

8.90 ±0.17 ±0.10 

US] 

15.36 

11.34 ±0.09 ±0.33 

HU 

53.00 

13.10 ±0.30 

HU 

4.721 

9.25 ±0.45 

m 

15.37 

11.20 ±0.60 

HU 

53.10 

13.00 ±0.30 

HU 

4.721 

9.16 ±0.37 


15.56 

11.30 ±0.20 ±0.20 

m 

62.30 

13.02 ±0.27 

HU 

4.934 

9.03 ±0.30 ±0.20 

m 

16.26 

11.60 ±0.40 ±0.17 

HU 

62.40 

13.30 ±0.30 

HU 

4.953 

10.80 ± 1.00 

0 

16.83 

12.20 ±0.50 

HU 

62.50 

12.80 ±0.20 

HU 

4.953 

9.24 ±0.35 

m 

16.83 

11.57 ±0.23 ±0.20 

m 

63.00 

13.30 ±0.60 

HU 

5.150 

10.32 ±0.17 ±0.30 

m 

16.83 

12.10 ±0.30 

HU 

200.0 

16.30 ±1.60 ±0.90 

HU 

5.562 

10.31 ±0.15 ±0.30 


16.91 

11.46 ±0.06 ±0.33 

HU 

7000 

19.89 ±0.02 ±0.27 

HU 

6.105 

10.24 ±0.11 ±0.30 


18.14 

11.63 ±0.09 ±0.33 


7000 

19.73 ±0.14 ±0.26 

HU 

6.171 

10.00 ±0.20 

m 

18.17 

11.80 ±0.50 

HU 

8000 

19.90 ±0.30 

HU 

6.522 

10.47 ±0.14 ±0.30 

m 

18.17 

11.52 ±0.11 ±0.20 

cu 




6.920 

10.48 ±0.13 ±0.30 


19.37 

11.56 ±0.12 ±0.20 






*The point is excluded from the fit procedure. 
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Table 9. Experimental Bpp for intermediate |t| domain 


yi, GeV 

Bpp, GeV-2 

Ref. 

yi, GeV 

Bpp, GeV-2 

Ref. 

yi, GeV 

Bpp, GeV-2 

Ref. 




Linear parameterization \n{da/dt) oc ( 

-B\t\) 



2.034 

1.60 ±0.60 


3.942 

7.37 ±0.29 


23.00 

10.30 ±0.20 

m 

2.037 

0.20 ± 1.40 


4.220 

7.75 ±0.11 

m 

23.76 

10.60 ±0.30 

EU 

2.058 

0.59 ±0.11 


4.331 

8.39 ± 20 

El 

30.81 

10.91 ±0.22 

m 

2.066 

1.00 ±0.50 


4.519 

8.95 ± 0.06 

m 

31.00 

10.92 ±0.15 

m 

2.075 

0.00 ± 0.70 


4.621 

8.11 ±0.17 ±0.10 

[S] 

31.39 

10.93 ±0.20 

m 

2.098 

0.90 ± 0.40 


4.694 

8.53 ± 0.22 

El 

44.61 

11.10 ±0.15 

EU 

2.112 

0.40 ± 0.60 


5.010* 

8.19 ±0.13 

m 

44.87 

10.83 ±0.20 

m 

2.132 

1.50 ±0.60 


5.122 

8.83 ± 20 

El 

45.00 

10.75 ±0.75 

EU 

2.142 

1.20 ±0.40 


5.780 

8.86 ± 23 


45.41 

11.08 ±0.20 

El 

2.153 

1.30 ±0.30 


6.028 

8.58 ± 0.24 

m 

46.34 

10.84 ±0.20 


2.166 

3.00 ± 0.50 


6.223 

9.12 ±0.19 

El 

52.61 

10.90 ±0.15 

EU 

2.180 

2.70 ± 0.50 


9.081 

9.65 ±0.07 

m 

52.80 

10.70 ±0.20 

EU 

2.191 

3.60 ± 0.60 


9.081 

10.20 ± 20 

m 

52.80 

10.34 ±0.19 ±0.06 

EU 

2 .201* 

5.10 ±0.60 


9.778 

9.40 ±0.10 ±0.14 


52.99 

10.80 ±0.20 

EU 

2.688 

5.55 ±0.20 

m 

10.69 

11.10 ±0.90 

m 

53.00 

11.06 ±0.11 

EU 

2.780 

6.46 ±0.25 


11.54 

9.50±0.10± 0.14 

HU 

53.00 

10.20 ±0.10 

EU 

2.837 

5.55 ±0.23 


16.26 

9.60 ±0.10 ±0.14 


53.35 

10.78 ±0.23 

El 

3.034 

6.25 ±0.26 


16.66 

10.38 ±0.17 


54.44 

10.93 ±0.18 


3.363 

7.12 ±0.11 

m 

16.66 

9.20 ± 0.50 

m 

62.00 

10.30 ±0.30 

m 

3.378 

7.37 ±0.25 

m 

19.66 

8.90 ± 0.70 

m\ 

62.00 

10.71 ±0.08 

EU 

3.433 

7.60 ± 0.50 

E] 

19.66 

9.90 ± 0.40 


7000 

23.60 ± 0.50 ± 0.40 

EU 

3.778 

7.94 ± 0.26 

m 

21.49 

10.42 ±0.17 

m\ 




3.916 

7.95 ±0.18 

m 

21.84 

10.38 ±0.22 

El 






Quadratic parameterization ln{da/dt' 

oc (- 

B\t\±Ct‘^) 



2.212 

4.60 ± 0.90 

E] 

3.826 

9.78 ±0.21 

m 

6.272 

9.17±0.11 

m 

2.236 

5.50 ± 0.70 


4.131 

7.00 ± 1.10 

El 

6.434 

9.50 ±0.90 

E] 

2.265 

6.40 ± 0.60 


4.220 

8.16 ±0.28 

m 

6.434 

9.80 ±0.30 


2.272 

7.64 ± 0.46 


4.220 

8.35 ± 0.25 


6.477 

10.90 ± 0.60 


2.304 

6.50 ±0.60 


4.286 

9.62 ± 0.22 

EU 

6.547 

9.63 ±0.78 

EU 

2.307 

8.12 ±0.84 


4.701 

9.79 ± 0.23 


6.871 

9.90 ± 1.00 

El 

2.311 

6.50 ±0.80 


4.729 

8.56 ± 0.47 

m 

6.929 

7.97 ± 1.56 

EU 

2.325 

7.70 ± 0.30 


5.010 

9.05 ± 0.34 

m 

7.138 

11.60 ±0.70 

El 

2.339 

6.30 ±0.50 


5.010 

9.71 ±0.16 


7.422 

9.10 ± 1.10 


2.342 

8.80 ± 1.00 


5.084 

10.03 ±0.28 

EU 

7.584 

9.34 ±0.28 

m 

2.360 

5.70 ± 0.50 


5.120 

9.90 ± 1.10 

El 

9.081* 

11.20 ±0.22 

m 

2.395 

6.30 ± 0.40 


5.128 

11.40 ± 1.20 


9.778 

10.30 ±0.10 ±0.15 

[lU 

2.458 

7.10 ±0.40 


5.440 

10.37 ±0.33 

EU 

11.54 

10.60 ±0.20 ±0.16 


2.510 

7.30 ± 1.00 


5.462 

8.89 ± 0.52 

EU 

13.76 

10.70 ±0.20 ±0.16 


2.682 

8.10 ±0.40 


5.491 

8.81 ±0.25 

m 

13.90 

11.40 ±0.70 

EU 

2.768 

6.10 ±0.90 


5.559 

10.00 ±0.80 

El 

16.26 

11.30 ±0.10 ±0.17 

HU 

2.768 

7.20 ±0.20 


5.757 

9.79 ± 0.40 

EU 

18.17 

11.30 ±0.10 ±0.17 


2.768 

7.80 ±0.15 

EH 

5.981 

9.60 ± 0.90 

El 

19.46* 

12.64 ±0.12 

HU 

2.972 

8.29 ±0.16 


6.028 

9.79 ± 0.63 

m 

27.60 

10.70 ±0.90 

EU 

2.978 

6.60 ± 0.70 

E] 

6.028 

9.96 ±0.21 


31.00* 

12.83 ±0.34 

EU 

3.077 

8.20 ± 1.00 


6.059* 

12.50 ±0.70 

El 

53.00 

10.09 ±0.47 


3.363 

8.46 ±0.16 

EH 

6.123 

9.06 ± 0.22 


62.00 

10.12 ±2.90 


3.497 

8.70 ± 0.50 

m 

6.212 

10.48 ± 0.43 

EU 




3.627 

8.46 ±0.16 

EH 

6.248 

8.68 ± 0.79 

EU 





‘The point is excluded from the fit procedure. 
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Table 10. Experimental Bpp for low |t| domain 


•/i, GeV 

Bpp 

GeV-2 

Ref. 

yi, GeV 

Bpp, GeV-2 

Ref. 

yi, GeV 

Bpp, GeV- 

2 

Ref. 

1.885 

83.70 ± 

24.00 

m 

1.946 

26.60 ± 3.20 

m 

2.638 

13.10 ±0.40 



m 

1.887 

129.0 ± 

25.00 

m 

1.946 

23.30 ± 1.10 

m 

2.669 

13.10 ±0.30 



0 

1.889 

61.20 ± 

16.50 

[58] 

1.948 

22.00 ± 2.00 

[60] 

2.957 

14.10 ± 1.60 




1.891 

62.30 ± 

7.40 

m 

1.949 

19.12 ±0.83 

m 

2.987 

12.90 ±0.39 

± 

0.08 

m 

1.891 

52.40 ± 

3.80 

[58] 

1.950 

20.90 ±2.10 

[63] 

3.077 

13.50 ±0.90 



0 

1.891 

71.50 ± 

4.50 ± 7.00 

m 

1.951 

25.90 ± 1.70 

m 

3.098 

12.80 ±0.70 

± 

0.08 

m 

1.892 

38.80 ± 

1.30 

ISO] 

1.951 

21.00 ± 1.10 

m 

3.524 

12.50 ±0.29 

± 

0.08 


1.894 

39.20 ± 

2.00 

m 

1.951 

24.70 ± 0.60 

[60] 

3.555 

12.20 ±0.39 

± 

0.08 


1.896 

51.20 ± 

4.60 

m 

1.956 

24.80 ± 1.20 

m 

3.612 

12.60 ±0.29 

± 

0.08 


1.896 

40.50 ± 

0.90 

[60] 

1.957 

22.20 ± 1.10 

m 

3.627 

13.50 ±0.90 



EQ] 

1.896 

47.70 ± 

2.70 ± 5.00 

m 

1.957 

24.10 ±0.60 

[SO] 

3.686 

12.20 ±0.59 

± 

0.08 

m 

1.898 

32.30 ± 

1.70 

[5H] 

1.960 

23.40 ± 0.60 


4.108 

12.66 ±0.05 



0 

1.900 

39.00 ± 

2.00 

m 

1.960 

22.80 ± 0.94 

m 

4.108 

13.00 ±0.50 



m 

1.901 

45.10 ± 

4.10 


1.962 

24.90 ± 1.20 


4.351 

13.00 ±0.10 



m 

1.902 

35.70 ± 

0.80 

m 

1.962 

18.79 ±0.80 

m 

4.540 

11.80 ±3.00 



0 

1.906 

42.20 ± 

4.80 

m 

1.969 

18.23 ±0.80 


4.540 

12.20 ±0.10 



US] 

1.907 

36.50 ± 

0.80 

[60] 

1.973 

20.40 ± 2.20 

IS] 

4.621 

11.38 ±0.17 

± 

0.10 

m 

1.908 

32.10 ± 

3.50 

m 

1.976 

21.73 ±0.85 

m 

11.54 

12.49 ±0.18 

± 

0.33 

ES] 

1.908 

33.50 ± 

1.20 

[E] 

1.985 

15.90 ± 0.80 

[SO] 

13.76 

13.20 ± 1.20 



m 

1.909 

44.15 ± 

5.55 

m 

1.989 

15.00 ± 1.80 

m 

15.36 

12.46 ±0.14 

± 

0.33 

ES] 

1.911 

28.60 ± 

6.80 

m 

1.994 

12.00 ±2.00 

[SO] 

16.91 

12.52 ±0.16 

± 

0.33 


1.912 

33.60 ± 

1.70 

[60] 

1.996 

14.10 ±2.60 

m 

18.14 

13.03 ±0.19 

± 

0.33 


1.917 

27.40 ± 

1.70 

m 

2.005 

16.30 ± 1.20 


19.42 

17.00 ±3.60 



m 

1.918 

24.80 ± 

2.60 

m 

2.006 

18.00 ± 0.50 

ISS] 

19.46 

12.68 ±0.10 

± 

0.32 

ES] 

1.918 

28.80 ± 

1.20 

[H] 

2.036 

18.60 ± 1.50 

m 

24.30 

12.30 ±0.50 

± 

0.07 

m 

1.920 

29.10 ± 

0.50 

[HD] 

2.107 

15.20 ±0.30 

[S3] 

30.40 

12.70 ±0.50 



m 

1.921 

22.70 ± 

4.20 

m 

2.115 

16.10 ±0.40 

m 

30.70 

12.60 ±0.30 



EH] 

1.922 

24.70 ± 

1.70 

[HD] 

2.138 

15.00 ± 0.30 


31.00 

11.37 ±0.60 



m 

1.923 

22.63 ± 

1.30 

m 

2.140 

14.90 ± 0.60 

m 

52.60 

13.03 ±0.52 



m 

1.926 

23.80 ± 

2.30 

m 

2.159 

14.20 ±0.30 

[60] 

52.80 

13.60 ±2.20 



[30] 

1.926 

26.80 ± 

1.20 

[SI] 

2.184 

15.20 ±0.30 


52.80 

13.92 ±0.37 

± 

0.22 

m 

1.926 

26.20 ± 

0.50 

m 

2.193 

15.40 ± 0.30 


52.80 

13.36 ±0.53 



[69] 

1.927 

23.38 ± 

1.15 

m 

2.204 

14.50 ±0.30 


62.30 

13.47 ±0.52 



ES] 

1.931 

22.30 ± 

1.50 

m 

2.221 

14.20 ±0.20 


62.50 

13.10 ±0.60 



m 

1.931 

23.70 ± 

0.40 

[BD] 

2.223 

14.20 ±0.30 

m 

540.0 

17.20 ± 1.00 



m 

1.932 

22.52 ± 

1.05 

m 

2.227 

15.80 ±0.30 

[60] 

540.0 

17.60 ± 1.00 



m 

1.933 

23.80 ± 

1.20 

[SI] 

2.254* 

5.20 ±0.30 


540.0 

17.10 ±1.00 



m 

1.936 

21.80 ± 

0.60 

[SS] 

2.260 

16.50 ±0.40 


541.0 

15.50 ±0.20 



m 

1.938 

18.80 ± 

1.90 

m 

2.292 

15.10 ± 1.80 

m 

546.0 

15.30 ±0.30 



m 

1.938 

18.98 ± 

0.90 

m 

2.294 

12.40 ±0.50 

[SO] 

546.0 

15.20 ±0.20 




1.939 

26.40 ± 

1.14 

[58] 

2.296 

14.40 ± 0.40 

[S] 

546.0 

15.28 ±0.28 

± 

0.09 

ES] 

1.940 

25.20 ± 

1.10 

m 

2.304* 

5.90 ± 0.30 

m 

1800 

17.20 ± 1.30 



m 

1.942 

22.40 ± 

2.70 

m 

2.325 

13.30 ± 1.30 


1800 

16.30 ±0.50 



Ea 

1.943 

20.01 ± 

0.82 

m 

2.348 

13.30 ±0.20 

m 

1800 

16.37 ±0.45 



m 

1.944 

20.10 ± 

2.20 

m 

2.352 

13.20 ±0.30 

m 

1800 

16.99 ±0.47 



ES] 

1.945 

22.80 ± 

0.40 

[60] 

2.430 

15.60 ± 1.40 

0 







‘The point is excluded from the fit procedure. 
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Table 11. Experimental Bpp for intermediate |t| domain 


yi, GeV 

Bpp, GeV-2 

Ref. 

yi, GeV 

Bpp, GeV-2 

Ref. 

yi, GeV 

Bpp, GeV-2 

Ref. 

Linear parameterization \n{da/dt) oc {—B 

t\) 

2.156* 

17.80 ± 0.40 

[5] 

4.307 

12.83 ±0.21 

0 

31.00 

11.16 ±0.20 


2.177 

11.80 ± 1.00 


4.540 

11.80 ±2.90 

m 

52.80 

10.68 ±0.20 

[36] 

2.271 

11.20 ± 1.30 


4.621 

12.09 ±0.24 ±0.10 

[SI 

53.00 

11.50 ±0.15 


2.335 

16.50 ± 1.00 


4.896 

12.33 ±0.79 

m 

62.00 

11.12 ±0.15 


2.465 

12.90 ±0.70 


4.896 

12.30 ±0.80 

0 

540.0 

13.30 ± 1.50 

[89] 

2.558 

13.40 ± 1.10 


4.934 

12.66 ±0.29 

m 

540.0 

13.70 ± 0.30 

m 

2.702 

13.00 ±2.30 


4.934 

12.70 ±0.30 

0 

540.0 

13.70 ±0.20 ±0.20 

m 

2.768 

13.00 ±0.80 


5.627* 

8.78 ± 1.00 

m 

546.0 

13.40 ± 0.30 

m 

2.857 

13.10 ±0.70 


5.642* 

11.44 ±0.20 

[82] 

546.0 

13.60 ±0.80 


3.550 

12.10 ±0.40 


5.642* 

11.64 ±0.08 

0 

546.0 

14.20 ± 0.40 


3.550 

12.60 ±0.20 

m 

6.621 

13.10 ±0.80 

m 

546.0* 

15.35 ±0.18 ±0.06 

[75] 

3.851* 

15.60 ±0.40 

[5] 

7.006 

11.80 ±0.10 ±0.12 

[84] 

1020 

16.20 ±0.50 ±0.50 

m 

3.923 

13.15 ±0.47 

m 

7.864 

13.26 ±0.30 

[H5] 

1800 

16.30 ±0.30 

m 

3.923 

13.10 ±0.40 

0 

7.875 

12.70 ±0.80 

m 

1800 

16.40 ±0.12 


4.108 

12.57 ± 0.20 

m 

8.777 

11.30 ±0.20 ±0.11 

m 

1800 

16.98 ±0.24 ±0.07 

[75] 

4.108 

12.72 ±0.11 

0 

11.54 

10.10 ±1.00 

m 

1960 

16.86 ±0.10 ±0.20 

m 

4.307 

12.84 ±0.21 

m 

13.76 

11.40 ±0.60 

m 






Quadratic parameterization ln{da/dt) oc 

{-B\t 

±Ct^) 

2.768 

12.20 ±0.80 

m 

8.777 

12.20 ±0.70 ±0.12 

m 

31.00 

13.06 ±0.39 

m 

2.972 

12.10 ± 1.00 


9.778 

12.60 ±0.20 ±0.19 

PS] 

53.00* 

10.96 ± 0.44 


3.363 

11.40 ± 1.00 


11.54 

12.80 ±0.30 ±0.19 


62.00* 

11.30 ±0.45 


3.627 

12.40 ± 1.50 


13.76 

11.90 ±0.50± 0.18 


1800 

16.25 ±0.23 

m 

6.621 

13.10 ±0.80 

[83] 

16.26 

12.60 ±0.40 ±0.19 


1800 

16.87 ±0.57 


7.006 

12.80 ±0.40 ±0.13 

m 

18.17 

13.10 ±0.40 ±0.20 


1800 

20.65 ± 3.84 


7.864 

13.26 ±0.30 

m 

19.46* 

13.27 ± 0.24 

m 





‘The point is excluded from the fit procedure. 












